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Abstract. 

The two-dimensional Potts model can be studied either in terms of the original 
Q-component spins, or in the geometrical reformulation via Fortuin-Kasteleyn (FK) 
clusters. While the FK representation makes sense for arbitrary real values of Q by 
construction, it was only shown very recently that the spin representation can be 
promoted to the same level of generality. In this paper we show how to define the 
Potts model in terms of observables that simultaneously keep track of the spin and FK 
degrees of freedom. This is first done algebraically in terms of a transfer matrix that 
couples three different representations of a partition algebra. Using this, one can study 
correlation functions involving any given number of propagating spin clusters with 
prescribed colours, each of which contains any given number of distinct FK clusters. 
For < Q < 4 the corresponding critical exponents are all of the Kac form hr.s, with 
integer indices r, s that we determine exactly both in the bulk and in the boundary 
versions of the problem. In particular, we find that the set of points where an FK cluster 
touches the hull of its surrounding spin cluster has fractal dimension d2.i = 2 — 2/i2,i. If 
one constrains this set to points where the neighbouring spin cluster extends to infinity, 
we show that the dimension becomes ^ — 2 — 2/11,3. Our results are supported by 
extensive transfer matrix and Monte Carlo computations. 



PACS numbers: 64.60.De 64.60.F- 05.50+q 



1. Introduction 



The two-dimensional Q-state Potts model [T] is one of the most well-studied models 
in the realm of statistical physics. It can be addressed on the lattice using 
traditional approaches [2] (duality, series expansions, exact transformations) and 
algebraic techniques [3] (Temperley-Lieb algebra, quantum groups). In some cases it 
can be exactly solved using quantum inverse scattering methods |4j. When < Q < 4 it 
gives rise to a critical theory in the continuum limit, which can in turn be investigated 
by the methods of field theory [5] (Coulomb gas, conformal field theory) or probability 
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theory [B] (Stochastic Loewner Evolution). Despite of this well-equipped toolbox and 
sixty years of constant investigations, the Potts model remains at the forefront of current 
research, most recently in the context of logarithmic conformal field theory [TllH]. 
The partition function of the Potts model reads 

^ = E n ' (1) 

(T {ij)eE 

where K is the couphng between spins (jj = 1,2, ... ,Q along the edges E of some 
lattice C. The Kronecker delta function Sa^^crj equals 1 if cxj = aj, and otherwise. The 
precise choice of C does not affect universal critical properties, so for convenience we 
choose the square lattice. The transition line — which gives rise to a critical theory for 
< Q < 4 — is then given by the selfduality criterion e^ = 1 + y/Q. 

There are two obvious ways to rewrite the local Boltzmann weights: 

where we have used that S^-^^aj = or 1. The partition function ([T]) is obtained by 
expanding a product of such factors; each term in the expansion corresponds to choosing, 
for each (ij) G E, either the first or the second term in ([2]). 

Consider first the expansion in powers of w~^. Each factor of comes with the 
term (1 — ^o-j.o-j) that forces the spins Uj and aj to take different values; there is thus a 
piece of domain wall (DW) on the edge (ij)* dual to (ij). Making this choice for each 
(ij) G E defines a DW configuration on the dual lattice C*. This DW configuration 
defines a graph G (not necessarily connected), the faces of which are the clusters of 
aligned spins. Since we do not specify the colour of each of these clusters, a DW 
configuration has to be weighted by the chromatic polynomial Xg*{Q) of the dual graph 
G*. Initially Xg*{Q) is defined as the number of colourings of the vertices of the graph 
G*, using colours {1,2, ... ,Q}, with the constraint that neighbouring vertices have 
different colours. This is indeed a polynomial in Q for any G, and so can be evaluated 
for any real Q (but Xg*{Q) is integer only when Q is integer). The partition function 
(II]) can thus be written as a sum over all possible DW configurations 

^spin=»^^$^t.-'<=^^*^(^)xG.(g), (3) 
G 

where N is the number of spins, and length(G') denotes the total length of the domain 
walls. 

Alternatively we may consider the expansion in powers of v. The set of edges 
(ij) G E for which the term vS^^^a^ is chosen defines a graph H, the connected 
components of which are known as Fortuin-Kasteleyn (FK) clusters [9]. Since spins 
belonging to the same FK cluster are obviously aligned, we can perform the sum over 
{ai} in ([1]) to obtain 

Zfk = J2 Q»"i"^*^'^'^(^)i;«'^<^ses(//) _ (4) 

H 
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We stress that spins belonging to different FK clusters can still be aligned (with 
probability 1/Q), even if the two clusters are adjacent on the lattice. 

Despite of the fact that Zspm = Z^vi = Z exactly, the different geometrical degrees 
of freedom {G resp. H) underlying the spin and FK formulations give rise to different 
critical exponents in the continuum limit for < Q < 4. In both cases it is possible 
to build Z from a transfer matrix that acts on respectively the FK [10] and spin [TT] 
degrees of freedom. The transfer matrix is an essential ingredient in both analytical 
and numerical studies, and it gives access to the rich algebraic structures underlying the 
Potts model. In the FK case this structure is the Temperley-Lieb algebra [3], whereas 
the particular partition algebra that emerges in the spin case [11] still awaits a complete 
study. 

Critical exponents corresponding to FK clusters were first derived by Coulomb gas 
techniques (see e.g. [5]) and have subsequently been confirmed by the rigorous methods 
of SLE (see e.g. [6]). In particular the probability that i distinct FK clusters propagate 
between small neighbourhoods of two distant points defines a critical exponent that can 
be found in the Kac table of CFT as Hq i in the bulk case, and as /ii,i+2^ in the boundary 
case |I2]. The spin case turns out to be richer, since one needs in addition to specify the 
respective colours of the propagating spin clusters. A pair of adjacent clusters having 
different (resp. identical) colours give rise to a thin (resp. a thick) DW. The exponents 
corresponding to ii thin DW and £2 thick DW can be identified as h^^^f.^^2h the bulk 
case and as /ti+2(£i-£2),i+4^i the boundary case [11]. At present this result is based 
on (substantial) numerical evidence, supported by analytical arguments for the special 
case £1 = 0. 

The purpose of this paper is to study the joint properties of FK and spin clusters, by 
simultaneously keeping track of the geometrical degrees of freedom G and H appearing 
in ([3])-(jl]). The corresponding geometrical observables are defined precisely in section [2j 
We show that making certain physical assumptions — and insisting on recovering the 
results for FK [12] and spin [11] clusters as special cases — leads to a conjecture for the 
scaling dimensions of joint observables in the general case. This conjecture generalises 
all the preceeding results and gives access to new information on the interaction between 
FK and spin clusters. All the scaling dimensions fit into the Kac table h^^s with integer 
indices r and s. 

This result constitutes a further step in the programme [13] of classifying non- 
unitary boundary conditions in two-dimensional geometrical models. 

In section [3] we define a transfer matrix that contains complete information about 
these joint spin-FK observables. It formally acts on three coupled partition algebras. 
We study its spectrum numerically in section HJ finding strong support for our general 
conjecture for the scaling dimensions of geometrical observables. In section [5] we further 
corroborate the conjecture by studying the most relevant among the new exponents by 
Monte Carlo simulations. We discuss our findings in section [61 
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Figure 1. Schematic example of excitations considered in this paper, with three spin 
clusters carrying different coulours propagating along the imaginary time direction. 
The shaded regions represent propagating FK clusters living inside the spin clusters. 
Using the notations ([5]), this excitation corresponds to 1^2^3'^. 

2. Joint spin and FK observables 

We begin by reviewing tlie definition of observables in tlie spin representation of tlie 
Potts model [11]. For tlie purposes of this discussion it suffices to consider two-point 
correlation functions. We then ask how the probability, that a certain number i of spin 
clusters — each corresponding to a given value of the spin cij — connect a small (in units of 
the lattice spacing) neighbourhood A to another small neighbourhood B, decays when 
the distance x between A and B increases. It is convenient to specify the two-point 
function by just giving the i colour labels (Tj. By the Sq permutation symmetry only 
the relative colours matter. For instance, with i = 4 the labels 1123 and 2241 specify 
the same correlation function. 

Alternatively one may think of this question in terms of the DW that separate 
pairs of adjacent (when moving along the rim of A or B) spin clusters. Each DW can 
be either thin or thick, depending on whether the two adjacent clusters have different or 
identical colours. The reason for this nomenclature is that adjacent spin clusters with 
different colours can touch (in which case the corresponding DW has a width of one 
lattice spacing), whereas if the colours are identical they cannot (since then there would 
be no separating DW). In other words, the width of a thick DW is at least two lattice 
spacings. The results of [11] show that this distinction is important, since it survives in 
the continuum limit and gives rise to different critical exponents. 

An obvious but important consequence of ([2]) is that if two points belong to the 
same FK cluster, it implies that they also belong to the same spin cluster. The opposite 
implication is however not true. This means that FK clusters live "inside" spin clusters, 
and in particular a given spin cluster can contain several distinct FK clusters. This 
suggests introducing a more general two-point function where each propagating spin 
cluster contains a given non-negative number of propagating FK clusters. We describe 
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this by a label 

(5) 

where cr^ is the colour of fc'th spin cluster, and = 0, 1,2, . . . is the number of FK 
clusters living inside the /c'th spin cluster. An example of such an excitation is shown 
in Fig. [U 

2.1. Result for bulk and boundary critical exponents 

Two different geometries are of interest. In the bulk case, the neighbourhoods A and B 
are at arbitrary, but widely separated, locations in the infinite plane. This is conformally 
equivalent to an infinitely long cylinder — a strip with periodic boundary conditions — 
with A and B situated at the two extremities. In the boundary case, the geometry is 
that of the upper half plane, with A being at the origin and B far away from the real 
axis. It suffices to consider free boundary conditions on the real axis, since the work of 
pT| has already demonstrated that fixed or mixed boundary conditions on the positive 
real axis can be accommodated by fusing the operator that inserts the propagating 
(spin and FK) clusters with an appropriate boundary condition changing operator at 
the origin. 

The results for the critical exponents can be stated in terms of the Kac 
parametrisation of CFT 

where 2 < k < 4 parameterises the number of states Q G [0, 4] via 

Q = 4(cos^)'. (7) 

In the bulk case (geometry of the plane) the critical exponent corresponding to the two- 
point function defined above is denoted h{Q) = h{Q; a^Vg^ . . . o""'^), and the correlation 
function (or probability) of the prescribed event decays like P oc x^^'^'^'^^ for x ^ 1. 
Equivalently, on a long cylinder of size LxM with M ^ L, periodic boundary conditions 
in the L-direction, and A and B identified with the opposite ends of the cylinder, the 
decay is exponential: P oc Q-^'^i^'^/^) KQ) _ the boundary case (geometry of the half 
plane) the critical exponent for free boundary conditions is denoted h{Q). 
We can write the bulk and boundary exponents in a unified way as 

h{Q) = V,,, , (8) 

HQ) = hl+2q,,l+2q2 , (9) 

where the right-hand side refers to ([6]) and we have defined a "charge" {qi,q2)- The 
result of [12] is then that i propagating FK clusters correspond to 

(gi, q^) = (0, i) , for FK clusters. (10) 

In the case of spin clusters with ii thin DW and £2 thick DW the result of [llj is that 

(?i) I2) = (^1 — £2, 2£i) , for spin clusters. (11) 
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These results strongly suggest that the costs for creating an excitation consisting 
of several types of interfaces simply add up at the level of the charge (gi, 52)- This is a 
well-known situation in the Coulomb gas (CG) formalism. But note that whereas the 
result (fTOj) for FK clusters indeed admits a CG derivation [12j no such argument has yet 
been established for the spin cluster result (fTTll . Admitting the additivity of charges, 
we can nevertheless read off the following rules from flTTj) 



We now argue that two additional rules are needed for dealing with the insertion 
of FK clusters inside a given spin cluster. Indeed, when a single FK cluster is inserted 
(ttj = 1) we create a pair of interfaces separating the hull of the FK cluster from the 
hull of the surrounding spin cluster. The insertion of each subsequent FK cluster (when 
ctj > 1) will in addition create an interface separating the FK cluster from the one 
preceeding it. To match ( ITOl) we clearly need (gi, 52) = (0, 1) in the latter case. To deal 
with the former, we note that using our conventions, a configuration with just i FK 
clusters can either be labelled 1^ (if each FK cluster has the same colour) or 1^2^ ... 
(if each FK cluster has a different colour). In either case, additivity of charges and 
consistency with ffTOj) requires the correct rule to be 



The main result of this paper is the following conjecture: 

• Consider the excitation labelled as in (|5]). Compute the corresponding charge 
(gi,g2) by adding up the contributions from the ak, using rules ( fT2l) . and those 
of the ak, using rules ( fT3l) . Then the bulk and boundary critical exponents are 
given by 

2.2. Remarks on the first spin cluster 

The reader will have noticed that the nature (thin or thick) of the DW are determined 
by the relative colours of two spin clusters. It therefore remains to discuss carefully how 
to account for the first one of the spin clusters making up the excitation. 

In the bulk case, when only a single spin cluster is present (£ = 1), it will almost 
surely wrap around the periodic direction (i.e., separate the neighbourhoods A and B). 
Our results do not apply to this case. Indeed, it is well known that the critical exponent 
in this case is the magnetic exponent /ii/2,o (resp. /io.oi) for a spin cluster with ai = 
(resp. «! > 1) propagating FK clusters, cf. (|TOl) . If, on the other hand, we impose the 
additional constraint that the spin cluster must not wrap around the periodic direction, 
our results do apply. The constraint then amounts to the presence of one thick DW that 
separates the spin cluster from itself. 

In the boundary case, our results always apply, provided that the first spin cluster 
is counted as a thin DW. This makes intuitive sense, since the free boundary conditions 




for a thin DW , 
for a thick DW . 



(12) 




1,-1) for the first FK in a given spin cluster , 
1) for each subsequent FK in that cluster . 



(13) 
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permit the spin cluster to touch the boundary, just as in the case of two spins clusters 
with different colours. 

2. 3. Possible values of the exponents hr,s 

It is obvious that all excitations of the type ([5]) will lead to charges {qi,q2), where qi 
and q2 take integer (positive or negative) values. We now wish to clarify precisely which 
values (gi, ^2) G 1? can be obtained through the joint spin-FK observables. 

The spin part of the excitation is described by the numbers (^1,^2) of thin and 
thick DWs, as in ( ITTl) . We first discuss the boundary case, for which the allowed values 
are ii > 1 (since the first DW is thin) and £2 > 0. Using rules (fT2!) - ffT3l) one can then 
deduce the possible values of the charge {qi,q2)- 

• For qi > 1, we have the constraint q2 > 2gi. Indeed, the excitation with 
(^1,^2) = (q'i^O) and no FK clusters gives the charge (gi,2gi), while the one with 
(^1,^2) = {qi + 1,0) and a single FK cluster on one of the spin clusters gives the 
charge (gi,2gi + 1). Any higher value of q2 can be attained by placing additional 
FK clusters on the same spin cluster. 

• For even gi < 0, we have the constraint q2 > ^^y^- Indeed, the excitation with 
(^1,^2) = (1, — ^) and one FK cluster on each one of the spin clusters corresponds 
to q2 = ^^Y^, while any higher value can be obtained by adding further FK clusters. 

• For odd qi < 0, we have the constraint q2 > ^3^- Indeed, the excitation with 
(^1, ^2) = (1, ~^^) and one FK cluster on each of the spin clusters except the first 
one saturates the inequality on q2, whose value can be increased by the addition of 
more FK clusters. 

The discussion of the bulk case is similar, except that the allowed values of (^1,^2) 
are £1,^2 > with £1 7^ 1 (because of the periodic boundary conditions), and of course 
(^1; ^2) 7^ (0, 0) in order to have a non-trivial excitation. The leads only to a very minor 
modification of the constraints on the charge (gi, ^2) derived in the boundary case: For 
gi < we have now g2 > [^1 • 

In conclusion, the spin-FK observables described in this paper make it possible to 
produce (say, in the bulk case) all the Kac table exponents h^^s with integer indices (r, s) 
above or below both of the lines s = 2r and s = The two cones in between these 

lines are not accessible by the spin-FK observables. 

3. Transfer matrix construction 

We now describe a transfer matrix whose spectrum contains all of the excitations 
discussed in section [2J 

I That is, up to discretisation effects due to the fact that (r, s) are integers — see above for precise 
statements. 
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First recall that a partition "P of a set X is a set of nonempty subsets of X such 
that every element z e X is in exactly one of these subsets. The elements of V (i.e., the 
nonempty subset of X) are called blocks. A block containing precisely one element of 
X is called a singleton. A partition Va is said to be a refinement of Vb, and we write 
Va ^ Vb, provided that each block in Va is a subset of some block in Vb. This defines a 
partial order of the partitions of X. 

We shall need a few simple operators acting on V. We denote the identity operator 
by /. For i,j G X the join operator Jij acts as / if i and j belong to the same block; 
if not, it amalgamates the block containing i with the block containing j so as to form 
a single block. The detach operator Di acts by detaching i from its block, i.e., by 
transforming it into a new singleton block {i}. 

If we associate a Potts spin ai with each i G X we can finally define the indicator 
operator as Aij = 5a„<T, ■ I- 

In the representation theory of the Temperley-Lieb algebra, Di is usually defined 
in a slightly different manner. Namely, Di acts as Q ■ J if i is a singleton (i.e., it gives a 
Boltzmann weight Q)\ otherwise it transforms i into a singleton (with a Boltzmann 
weight 1). Then, setting e2i-i = Q^^^'^Di and e2i = Q^^^^i.i+i, the provide a 
representation of the Temperley-Lieb algebra [3] corresponding to the Q-state Potts 
model. In our construction of the transfer matrix we shall account for Q in a different 
manner, which is why we have set Q = 1 above. 

3.1. States and elementary transfer matrices 

The states of the Potts model in the spin-FK representation are given by a triplet of 
partitions {Vi,V2,V^) of the set X of vertices within a row of the lattice C 

The first partition Vi is arbitrary (not necessarily planar) and describes the spin 
colours: we have ai = aj if and only if i,j belong to the same block in Vi. The second 
partition V2 is planar and describes the connectivity of spin clusters. We have V2 :^Vi, 
since spins in the same spin cluster have equal colours. The third partition V3 is again 
planar and describes the connectivity of FK clusters. We have V3 :<V2, since spins in 
the same FK cluster are also in the same spin cluster. 

We represent graphically the transfer matrix as 

T = 

"' (14) 

The elementary transfer matrix Th that adds a horizontal edge between vertices i and j 
acts in Vi®V2® V3 as follows: 

Th = 1 ■ (/ - Aij) (g) / O / + (e^ - 1) ■ Aij (g) Jij Jij + 1 ■ A,j ® Jij ® / , (15) 

where the number before the dot (■) is the Boltzmann weight. The first term corresponds 
to the two spins ai and aj being different, in which case i,j can neither be in the same 
spin cluster, nor in the same FK cluster. The second term corresponds to i,j being 
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in the same FK cluster, hence also in the same spin cluster. Finally, the third term 
corresponds to i,j having the same spin without being in the same FK cluster; in this 
case the spin clusters must be joined up since i,j are neighbours on the lattice. 

We can similarly define the elementary transfer matrix that adds a vertical edge 
edge to the lattice, by propagating the vertex z to a new vertex i'. It acts in Vi®V2®Vz 



and each term has the same interpretation as in Th. Since T builds up the partition 
function Z, it must also account for the summation over the spin varibles, c/. ([1]). 
The easiest convention is to let the action of be accompanied by a sum over a^/. 
This sum is dealt with in the same way as in [TT]. Namely, let S denote the set of 
spin values {uj} being used in a given state, so that the number of elements |S| is just 
the number of blocks in Vi. After summing over dj/, the second and third terms in 
( IT6|) each correspond to a single non-zero contribution. The first term gives |S| — 1 
contributions where cTj' 7^ (Xj, but (Xj/ G S. The remaining contributions, where cTj' 7^ cTj 
and (jj/ ^ S can be regrouped (using the overall Sq symmetry) as a single contribution, 
with Boltzmann weight Q — in which {i'} becomes a singleton in V[. It is precisely 
because of this regrouping that it now makes sense to promote Q to an arbitrary real 
variable [TT]. 

The complete row-to- row transfer matrix T for a system of size L, shown 
graphically in f|T^ . is then obtained as the product of the elementary transfer matrices 
corresponding to all horizontal and vertical edges in a row: 



where L' = L — 1 in the boundary case (strip geometry) and L' = L with indices i 
considered modulo L in the bulk case (cylinder geometry). 
Several remarks are in order: 

(i) When representing Vi one can replace the actual spin values (jj by colour labels 
Cj defined such that q = Cj if and only if ai = aj. Using the Sq symmetry these 
colour labels can be brought into a standard form, thus reducing the number of 
basis states. Details of this construction are given in [IT) . 

(ii) For a system of size L, at most L different colour labels are used. Therefore the 
number of basis states depends only on L, and not on Q. In particular the state 
space is finite. 

(iii) The three partition algebras Vi, V2 and V3 are coupled by virtue of (fT5l) - (fT6|) . It 
is an interesting problem, beyond the scope of this paper, to analyse in details this 
situation from a representation theoretical point of view. 

(iv) If one neglects the information on FK clusters contained in V3, the second and third 
terms in (fT5!) - (fT6l) can be resummed, and one recovers the construction of [TT] . 



hke 



= 1 ■ (J - An,) ®Di®Di + (e^ - 1) ■ + A^. ®I®D,, (16) 




(17) 
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(v) If one neglects the information on spin clusters contained in V2, the first and third 
terms in (fT5|) -f lT6|l can be resummed. The result is then independent of Vi and 
reproduces the well-known Temperley-Lieb representation [3]. 

(vi) In the above construction we have accounted for Q in the partition algebra Vi. 
Therefore, the detach operators Di acting on V2 and have parameter Q = 1. 
In particular, the FK clusters described by V3 can be thought of as percolation 
{Q = 1) clusters inside the spin clusters described by Vi and V2- Leaving Q as a 
free parameter would amount to studying FK clusters of a second Q-state Potts 
model defined on top of the spin clusters of the original Q-state Potts model. This 
appears to be an interesting way of coupling a pair of (Q, Q)-state Potts models — 
and we intend to report further on this elsewhere. 

The leading eigenvalue Aq of the transfer matrix T gives the ground state free 
energy /q = — -^logAg. This /o coincides precisely with that of the usual FK transfer 
matrix ^0\, even when Q is non- integer. Its finite-size corrections possess a universal 

term whose coefficient determines the central charge of the corresponding CFT |14j . 

3.2. Correlation functions 

To obtain the two-point correlation defined in section [2] we need a variant transfer 
matrix T' which imposes the propagation of the defect labelled as in ([5]) along the 
(imaginary) time direction of the cylinder (or strip). From its leading eigenvalue Aq 
one can determine the energy gap A/ = —j^ log(AQ/Ao) whose finite-size scaling in turn 
determines the critical exponents h{Q) and h{Q) [141 [5]. 

To construct T' we modify the basis states Pi ® P2 ® V3 by marking some of the 
blocks in the partitions V2 and V3. First, we mark i blocks in V2 whose spin colours 
in Vi coincide with the choice of {ai, 0-2, ... , ai} in (E]). Second, let V3 be a refinement 
of V2 such that the fc'th marked block in V2 is refined into at least blocks in V3, of 
which precisely are marked. 

In order to conserve the marked clusters in the transfer matrix evolution, none of 
the marked clusters must be "left behind", and two distinct marked clusters must not 
be allowed to link up. Imposing these rules in a precise way is tantamount to defining a 
modified action of the join and detach operators, Jij and Di, on the marked basis states. 
Namely, Jij is modified to give a zero Boltzmann weight if i,j are in distinct marked 
blocks (this prevents a marked block from disappearing), and Di is modified to give a 
zero Boltzmann weight if {i} is a marked singleton block (this prevents marked blocks 
from being "left behind"). Finally, when Jij joins a marked and an unmarked block, 
the result is a marked block. 

The full state space corresponding to the excitation ([5]) is generated by letting T' act 
a sufficient number of times on a suitable initial basis state. This initial state is such that 
the fc'th marked block in V2 consists of consecutive points {ik, ik + 1, . . . ,ik + Oik — l} 
in a row of the lattice. In the refining partition V3 each of these points is a marked 
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1120 



122° 



Exact 



2 
3 
4 
5 



3.0000(1) 
5.005(5) 
6.97(2) 
8.85(10) 
2p-l 



5.01(1) 
8.05(5) 
11.02(3) 
13.9(1) 

3p- 1 



— 3.99(1) 

9.03(2) 7.03(3) 

11.2(1) 9.95(5) 

13.4(3) 12.7(3) 

2p + 3 3p-2 



Table 1. Bulk critical exponents corresponding to four different sector labels (O, as 
functions of the parameter p — The conjectured exponents read /ii^s for sector 

1^2°, hi^i for sector 1^2°, /i3,i for sector l^fO, and /i2,5 for sector l^2°?P . The table 
entries give the value of \p\, when (|6|) is rewritten as hr^s = {p^ — l)/(4p(p + 1)), with 
error bars shown in parentheses. 



singleton. Finally, those of the L points which are not marked in this construction of 
the initial state are taken as singletons, both in V2 and in P3. 

In summary, the modified transfer matrix T' keeps enough information, both about 
the mutual colouring of the sites and about the connectivity of spin and FK clusters, to 
give the correct Boltzmann weights to the different configurations, even for non-integer 
Q, and to follow the time evolution of the excitation defined by ([5]). 

4. Exact diagonalisation results 

We have numerically diagonalised the transfer matrix in the spin-DW representation for 
cylinders (resp. strips) of width up to L = 8 (resp. L = 9) spins. We verified that the 
leading eigenvalue Aq in the ground state sector coincides with that of the FK transfer 
matrix, including for non-integer Q. As to the excitations Aq, we explored systematically 
all possible colouring combinations (|5]) for up to ^ = 3 marked spin clusters with different 
choices for the number of FK clusters ak-, for a variety of values of the parameter k. 

Finite-size approximations of the critical exponents h{L) and h{L) were extracted 
from the leading eigenvalue in each sector, using standard CFT results [HI [5], and 
fitting both for the universal corrections in and the non-universal term. These 
approximations were further extrapolated to the L — 00 limit by fitting them to first 
and second order polynomials in gradually excluding data points corresponding to 
the smallest L. Error bars were obtained by carefully comparing the consistency of the 
various extrapolations. 

Representative final results for the bulk case (periodic boundary conditions) are 
shown in Table [H Corresponding results for the boundary case (free boundary 
conditions) are given in Table El In all cases the agreement with the conjecture made 
in section |2] is very good. 
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112° 




1^0 


II2030 


2 


5.001(2) 


9.00(1) 


9.01(2) 


7.002(4) 


3 


9.01(1) 


15.0(1) 


13.02(2) 


13.05(8) 


4 


12.92(5) 


20.98(2) 


17.02(3) 


18.85(10) 


5 


16.7(2) 


26.7(2) 


21.05(10) 


24.5(3) 


Exact 


4p-3 


6p — 3 


Ap+1 


6p — 5 



Table 2. Boundary critical exponents with free boundary conditions, corresponding to 
four different sector labels ([5]), as functions of the parameter p = The conjectured 
exponents read ft-sj for sector 1^2°, ^13^9 for sector 1^2", /i_i,3 for sector 1^1'^, and 
/i5_ii for sector 1^2*^3°. The table entries give the value of \p\, when ([6]) is rewritten as 
hr,s — {p^ — l)/(4p(p + 1)), with error bars shown in parentheses. 




(a) (b) (c) 

Figure 2. Various critical curves in the 3-state Potts model. The three pictures 
show the same configuration with different geometrical features emphasised, (a) Spin 
cluster configuration. There is a non-trivial (i.e. non-homotopic to a point) red cluster 
wrapping along one periodic boundary condition, but not the other one. One of its 
boundaries is shown in red, and contributes to the fractal dimension dgpin — 2 — 2hifl. 
(b) Same spin configuration, but this time the FK clusters and their corresponding 
boundaries are shown. The non-trivial red cluster is broken into smaller FK clusters, 
one of them remaining non-trivial. One of the two boundaries of the percolating FK 
cluster is shown in yellow and has a fractal dimension dpK = 2 — 2/io,i in the continuum 
limit, (c) Superposition of the previous red and yellow interfaces. Their intersection, 
shown in orange, corresponds to the set of points where the FK cluster touches the 
hull of its surrounding spin cluster. Its fractal dimension is (i2,i = 2 — 2ft,2,i- 

5. Monte Carlo simulations 

Although we the result of section [2] apply to a general excitation (|5]), only a few of 
the resulting scaling dimensions hr^s are relevant (i.e., < hr^s < 1) for some Q in the 
interval < Q < 4. Moreover, some of the excitations are relevant only in a part of the 
interval, where there are not "enough colours" to realise the choice of different cjfc used 
in (jS]). An example of this situation is the bulk operator 1^2°3° with scaling dimension 
/i2,5, which is relevant only for < Q < 2, meaning that Q is not large enough to 
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accommodate the three different spin colours defining the excitation. 

These remarks become important if we want to measure the fractal dimension 
corresponding to a given excitation in a Monte Carlo simulation. For convenience we 
restrict the discussion to bulk critical exponents. Among the excitations which involve 
both spin and FK degrees of freedom, it appears that only two are relevant and physical 
(in the above sense). The first of these is 1^ with wrapping around the periodic direction 
disallowed. It describes the insertion of a spin cluster with one FK cluster inside it. The 
charge is {qi,q2) = (—1)0) + (—1, —1) from ( |T2l) -( |T3|) leading to the scaling dimension 
/i_2,-i = ^2,1 by application of ([8]). The corresponding codimension 

d2,i = 2 - 2/i2,i = 3 - - (18) 

K 

is thus the fractal dimension of the set of points where the FK cluster touches the hull 
of its surrounding spin cluster. We have < i < 2 for all < Q < 4. An example of 



a geometrical curve corresponding to the dimension d2,i is shown in Fig. 2(c) 

The second excitation of interest is 1^2°. The charge is now {qi,q2) = 2 x (1,2) + 
(—1,-1) from (fT2l) -(!T3ll. and the corresponding scaling dimension reads hi^s by (jS]). 
The fractal dimension 

rfi,3 = 2 - 2/ii,3 = 4 - K (19) 

describes the same set of points as above, but with the additional constraint that a 
second spin cluster (the 2° part of the label) adjacent to the one whose hull is being 
touched by its internal FK cluster (the 1^ part of the label) must now also propagate all 
the way to infinity. Note that di^s and d2^i should coincide for the Ising model {Q = 2) 
and indeed we find di^^ = d2^i = 1 in this case. 

We checked numerically these fractal dimensions using Monte Carlo (MC) methods. 
Efficient MC algorithms for the non-integer Q-state Potts model are already available. 
We chose here to work with the Chayes-Machta algorithm [15] that works for Q G [1, 4]. 
This algorithm is easy to implement and allows one to keep track of both FK and spin 
clusters, even for non-integer Q [16]. 

The elementary step of the algorithm reads: 

• Find all the FK clusters in the configuration. 

• Independently label the clusters as active or inactive with respective probability ^ 
and Sites belonging to an active clusters are said to be active. 

• Erase all bonds. Independently add bonds between pairs of active sites with 
probability pc = I — = ■ The resulting clusters are the new FK clusters. 

Active sites correspond to spins in a given colour Qq, so that the spin clusters of this 
colour can be obtained by performing the bond-adding step of Chayes-Machta at zero 
temperature, i.e. by replacing pchj Pq = 1. The clusters and their boundaries are then 
detected in a standard way. 

Since the fractal dimensions we wish to measure are bulk properties, we work with 
an L X L lattice with periodic boundary conditions. We perform the statistics on 
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Ql : : : : : : : , 

0.5 1 1.5 2 2.5 3 3.5 4 



Q 

Figure 3. Fractal dimensions dpx, rfspin, ^2,1 and di^s (see text for details) as functions 
of Q. The small circles with their error bars are obtained from Monte Carlo simulations, 
while solid lines represent the CFT prediction. 

non-trivial clusters wrapping once around one of the periodic boundaries. The fractal 
dimensions dp are then obtained by measuring the mean value i of the curve length as 
a function of L = 32, 64, . . . , 800; using the relation i ~ AL'^^ . 

This algorithm allowed us to recover as a check the well-known fractal dimensions 
of the interfaces of the FK and spin clusters, which read respectively dpK = 2 — 2/io,i 
and (ispin = 2 — 2/ii q. The dimension d2,i = 2 — 2/i2,i of the set of points where the 
FK cluster touches the hull of its surrounding spin cluster can be measured without 
much more complication for Q > 1. To evaluate the dimension di^s = 2 — 2/11.3, we 
restrict the set of points considered when measuring 6/2,1 to the points whose adjacent 
spin cluster also wraps around the periodic boundary condition. As the Chayes-Machta 
for non-integer Q only keeps track of one spin cluster, we were able to measure c/1,3 only 
for Q integer. Note also that c/1,3 makes sense physically only for Q >2. 

Results are shown in Fig. [3l The agreement with our predictions is very good. The 
small deviations we observe when Q becomes close to Q = 4 are expected, as logarithmic 
corrections are known to occur in this case. 

6. Conclusion 

We have defined geometrical observables that keep track of both FK and spin clusters 
for any real Q G [0,4]. We conjecture that such observables are conformally invariant 
and we provide exact formulas for the bulk and boundary critical exponents. Our results 
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are supported by extensive transfer matrix and Monte Carlo computations. It is quite 
remarkable that all the exponents are of the Kac form hr^, an analytical understanding 
of those results would probably shed some light on this rather intriging point. We note 
also that we now have enough observables to cover all the Kac table h^^s for any integer 
choice of (r, s), except in the two cones delimited by the straight lines s = 2r and s = r/2 
(see section [2T3|1 . This remark is particulary important in the context of Logarithmic 
CFT (LCFT), where including more involved geometrical observables in the theory 
might yield some possibly unknown, interesting logarithmic features. In particular, our 
transfer matrix T at logarithmic points should have a much more complicated structure 
that the ones considered so far in lattice regularisations of LCFTs (see e.g. [3 E]). 
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